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We construct a class of solutions to the Einstein's equations for dimensions greater than or equal 
to six. These solutions are characterized by a non-trivial warp factor and possess a non-compact 
extra dimension. We study in detail a simple model in six dimensions containing two four branes. 
One of each brane's four spatial directions is compactified. The hierarchy problem is resolved by 
the enormous difference between the warp factors at the positions of the two branes, with the 
standard model fields living on the brane with small warp factor. Both branes can have positive 
tensions. Their positions, and the size of the compact dimension are determined in terms of the 
fundamental parameters of the theory by a combination of two independent and comparable effects — 
an anisotropic contribution to the stress tensor of each brane from quantum fields living on it and a 
, contribution to the stress tensor from a bulk scalar field. One overall fine tuning of the parameters 

of the theory is required — that for the cosmological constant. 

s 
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■ i. introduction 
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The work of Arkani-Hamed, Dimopoulos and Dvali (ADD) JlJ and of Randall and Sundrum (RS) 0] has stimulated 
interest in explaining the observed weakness of gravity (the "hierarchy problem") using extra dimensions. The ADD 
solution requires the new dimensions to have finite but large volume, which introduces a new hierarchy between the 
, volume of the compact dimension and the fundamental scale of the theory^. RS proposed instead five dimensional 
spacetime with curvature comparable to the fundamental scale, and showed that a massless graviton can be localized 
0^ ! to a 3+1 dimensional hypersurface known as the "Planck Brane". In their setup the extra dimension may be taken 
0^ ■ infinitely large, and four dimensional general relativity still agrees to high precision with long distance experimental 
"p^ ' measurement. The weakness of observed gravity is explained provided the standard model fields are localized to a 3+1 
dimensional "TeV Brane" , where the graviton wave function is small. Due to the exponential fall-off of the graviton 
O ' wave function away from the Planck brane, the distance between the TeV brane and the Planck brane does not need 
to be large in units of the fundamental scale. Randall and Lykken (RL) |^] have shown that an infinite fifth dimension 
i-£-h i is experimentally quite consistent with such a resolution of the hierarchy problem. 

In such a picture, it is necessary to introduce dynamics which determines the location of the TeV brane relative to 
the Planck brane. If this interbrane distance is not fixed, it becomes a massless modulus which leads to unnacceptable 
cosmology || and experimental consequences. Goldberger and Wise showed that adding a scalar field which propagates 
in the bulk and has a source on the branes is sufficient to fix this distance [Q. Several other suggestions have been 
made for bulk dynamics to fix the extra dimensional configuration j|] . 

In this paper we first construct a class of solutions to the Einstein's equations for dimensions greater than or equal 
to six. These solutions are characterized by a non-trivial warp factor and possess a non-compact extra dimension. 
We then study in detail a simple model in six dimensions containing two four branes that employs a metric of this 
form to address the hierarchy problem. One of each brane's four spatial directions is compactified on a circle of small 
radius. The hierarchy problem is resolved as in the RS and RL models by the enormous difference in the warp factors 
at the positions of the two branes, with the standard model fields living on the brane with small warp factor. 

In this model the positions of the branes, and hence the magnitude of the hierarchy, are determined by the combi- 
nation of two independent effects. The first is an anisotropic contribution to the stress tensor of each brane arising 
from the quantum effects of fields localised to it. The second is a contribution to the stress tensor from a scalar bulk 
field as in the model of Goldberger and Wise. These effects are naturally comparable in size and together can yield a 
sufficiently large value of the brane spacing to solve the hierarchy problem without fine tuning of parameters. 
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1 For examples of theories which naturally have finite but exponentially large volume for the additional dimensions see ref. Q. 
For earlier work on large extra dimensions and/or a low quantum gravity scale see refs. W. 



1 



The reason for the first effect is that the theory contains a compact dimension, in addition to the noncompact 
dimension r. The size of the compact dimension is in general an r dependent function, which is determined from the 
Einstein's equations. We argue that in general the component of the brane tension in the compact dimension will 
depend on its size, due to the quantum effects of fields localized to the brane. Then a consistent solution to Einstein's 
equations will fix each brane location at a particular value of r. However in this model this effect by itself does not 
give rise to a large hierarchy without fine tuning. Nevertheless when combined with the effect of a bulk scalar field 
on the geometry it is possible to realise a large hierarchy without fine tuning of parameters. 

Both branes in the theory can have positive tensions and the solution is free of singularities where general relativity 
might break down. One overall fine tuning of the parameters in the theory is required to adjust the four dimensional 
cosmological constant to zero. 



II. THE ANSATZ, EQUATIONS OF MOTION AND SOLUTIONS 



We begin by looking for solutions to the Einstein equations in D dimensions, where D is greater than or equal 
to six, in the presence of a constant background bulk cosmological constant. We assume all sources other than the 
bulk cosmological constant are restricted to subspaces of lower dimension. Hence our approach will be to first solve 
the equations of motion in the bulk to obtain solutions with a number of constants of integration that can then be 
adjusted to find solutions for various boundary conditions. 

The action in the bulk is 

S = J d D xy/^G(2M°- 2 R - A B ) . (i) 

We label a general coordinate by x M where M runs from to (D-l). We restrict our search to metrics of the simple 
form 



ds 2 = f{z)r lllv dx^dx v + s(z)dy 2 + dz 2 (2) 

where fi and v run from to D-3. The remaining two coordinates are labelled by y and z. Here the warp factors / 
and s are assumed to be functions only of z. 

The Einstein's equations in the bulk take the form 

2M?- 2 (Rmn ~ 9mnR) = - ^9mnAb ■ (3) 
The nontrivial components of this equation are 



V,, m 1/ , f ,s' (3~J>) , 1,/VV f ff'Y f D 2 -9D + 18 \ 2 (D-2) 2 + (D-2) 
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where a is defined by 
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To solve these equations note that we can rewrite eqn. (^|) in the form 
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This has the form of a first order differential equation for y . This differential equation is straightforward to solve, 
and we can then obtain / itself by performing a simple integration. We then use the result obtained for / in eqn. (|^) 
and the problem of determining s also then reduces to performing a simple integral. The results are 
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Here /o, so and c are constants to be determined by boundary conditions. In the limits of vanishing c and infinite 
c we recover the usual anti deSitter (AdS) metric. It is easy to see that in fact for any values of these constants the 
warp factors / and s change very rapidly as function of z. In particular there are always values of z where they are 
changing exponentially quickly. This suggests that these metrics are good candidates for a solution to the hierarchy 
problem. 

It is possible to use the class of metrics above to find solutions to the Einstein equations for various source configu- 
rations and geometries^. In the next section we exhibit a potential solution to the hierarchy problem based on these 
metrics. 

III. THE MODEL 
A. The Metric 

In this section we limit our interest to a solution where the fifth dimension y is compact and corresponds to an angle 
in the higher dimensional space. We relabel y by <j) in this section and hereafter to emphasize its angular character. 
The angle <j) runs from zero to 2ir. We allow the coordinate z to be non-compact and run from zero to infinity, it 
corresponds to a 'radius' in the higher dimensional space. We relabel it by r to emphasize its radial character. 

The geometry of our model consists of two four branes localized in the higher dimensional space at different values 
of r. Their positions are specified by the equations r = a and r = b where a < b. They can therefore be thought of 
as being similar to the surfaces of two infinitely long 'concentric cylinders' in the higher dimensional space, with the 
regular four dimensions parallel to the common axis of the cylinders, the fifth dimension going around the surface but 
perpendicular to the axis, and the sixth dimension being the radius. (This intuitive picture does not account for the 
fact that the space is curved.) The standard model fields are localized to the brane at r — a. The hierarchy problem 
will be resolved by the enormous difference between the values of the warp factor at r = a and r = b. 

Because the coordinate <j> is compact the four brane appears as a three brane at sufficiently long length scales. 

The branes divide the space into three distinct sections; 0<r<a,a<r<b, and r > b. In general there is no 
reason for the bulk cosmological constants in these three sections to be the same since the branes may be separating 
different phases of the theory. In what follows we will assume that they are different and will associate the three 
regions with three different values of a; a%, cti and respectively. 

The solutions of the Einstein equations in the three regions will be of the form of equations (^|) and (|l^) but with 
different values of the constants /o, sq and c. We will give these constants an additional subscript i, where i runs 
from 1 to 3, in order to differentiate them in the three different regions. 

Now the constant C\ is fixed to be -1 by the requirement that the solution be non-singular at the origin. This 
requirement also fixes the value of Sol to be (2 16 / 5 /(25a 2 )). To see that this choice does indeed smooth out the 
singularity at the origin we first examine the behavior of the the functions f(r) and s(r) as r tends to 0, 

f(r) = const +0(r 2 ) (11) 
s(r) = r 2 + 0{r 4 ) . (12) 

We then go to the 'cartesian' coordinate system which is smooth at the origin. 

x = r cos (13) 
y = r sin (f> . (14) 



"Note that these solutions are coordinate transformations of the bulk solutions found in ref. 
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It is straightforward to verify that in this coordinate system the components of the metric and their first and second 
derivatives (which go into the Ricci tensor) are smooth at the origin, showing that there is no singularity there. 

The requirement that the metric be bounded at infinity determines that the space outside r = b be AdS. This 
corresponds to setting C3 to infinity while keeping the products C3/03 and c 3 s 3 finite. We relabel these products by 
/ 3 and s 3 respectively. 

We are now in a position to write down the forms of the solutions in the three regions. For r < a, 



/ = /oie Qir [l + e- 



(15) 



s = s 01 e Qir [l + e -* air ]" 5 " [1 - e-s ai T • 



(16) 



For b > r > a, 



f = fv2e a2r [l-c 2 e-^ r Y : 



s = s 02 e Q2r [l - c 2 e-2^]-[l + c 2 e 



For r > b, 



f = he 



(17) 
(18) 

(19) 



s = s 3 e 



(20) 



The constant /01 is determined by normalizing the warp factor / to 1 at the position of the visible brane. The 
constant S01 was determined earlier. All the other constants above as well as the positions of the branes a and b must 
be determined by matching across the branes. 

We now write down the action for the branes. 



S b 



-I 



dx 



A I 



-G(6(r - a)[L M i - Ai] + 6{r - b))[L M2 - A 2 ]) . 



(21) 



Here G is the determinant of the metric tensor in the five dimensional subspace, Ai and A 2 are the cosmological 
constants on the two branes, and L M \ and L M2 are the Lagrangians of matter fields localized to the branes. 
The stress tensor for each brane has the form 



Tai 



AB 



+ < T M AB > 



(22) 



where T A is the contribution from the cosmological constant and < T M > the expectation value of the stress tensor 
of the matter fields living on the brane. We will be working in the semi-classical limit, treating gravity completely 
classically but accounting for the quantum effects of matter localized to the branes. 

Assuming the matter on the brane is in its ground state, Tab is constrained by four dimensional Lorentz invariance 
to be of the form 
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(23) 



where f3 2 and 7 2 are constants associated with each brane that we require to be positive. The form of T A is more 
constrained since it is proportional to Gab , which would imply that 1 = j 2 if the only contribution to T came from 
the cosmological constant on the brane. Thus the deviation of T from the Gab form is due entirely to the contribution 
from the matter Lagrangian. In a subsequent section we will show that the contribution to the stress tensor from the 
zero point energies of fields living on the brane lead to (i ^ 7. In general (3 and 7 will depend on the geometrical 
factors a, b and c 2 due to the quantum contribution to the stress tensor from matter localized on the brane. 

Now the Einstein equations for the upper 5 by 5 block of the Ricci tensor get modified in the presence of the branes 

to 
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2M*(R AB - g AB R) = -~g AB A B + 8nS{r - a)T% B + 8^(r - b)T b AB . (24) 

Since we already have the solution in the bulk we can get the complete solution by matching across the brancs. 
The metric tensor is continuous across the branes but its derivatives are not. The above equation fixes the jump 
discontinuity in the derivatives across the boundary. 

In terms of components the conditions on the derivatives are 

_^A^_I A -=/? 2 (25) 
2 / 2 s H y ' 

2A^r = - 7 2 . (26) 

We wish to apply these conditions to our solution. To simplify matters we first define 

f = A0 2 - 3 7 2 (27) 
9n>r 

F(c,r,a) = T - 5 v . (28) 

I e2™ — c) 

The conditions on the derivatives at the first boundary are 

1 



2 
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a 2 - ot\ + F(c 2 ,a, a 2 ) - F(-l, a, a{) = - ~7 X (29) 



a 2 - oti - - [F(c 2l a, a 2 ) - F(-l, a, on)] + - [F(-c 2 ,a, a 2 ) - F(l,a, ax)] = --jli ■ (30) 
The conditions at the second boundary are 

- a 3 - a 2 - F(c 2 ,b, a 2 ) = --7I (31) 

3 5 1 

- a 3 - a 2 + ^F{c2, b, a 2 ) - -F(-c 2 , b, a 2 ) = - -72 2 . (32) 

Looking at the above equations we have three parameters c 2l a and b which have to satisfy four independent 
equations. Hence a fine tuning is necessary. This is the fine tuning necessary to set the effective four dimensional 
cosmological constant to zero. Once this fine tuning has been made it is straightforward to find solutions to the above 
equations without singularities where all parameters are of order one in terms of the fundamental scale M* . However 
since we want to generate a hierarchy we want a 2 b — O(40) 3> 1. From now on we will assume that b is the only large 
parameter in the problem and that it is therefore responsible for generating the hierarchy. Subtracting eqn. ( pl| ) from 
eqn. (|3|) we see that 

^F(c 2 , b, a 2 ) - |F(-C2, b, a 2 ) = -\p 2 2 + ]pl . (33) 

Note that from the definition of 7 it is clear that if p = 7, then 7 = 7. But the difference between /3 and 7 arises 
from the vacuum energy of quantum fields localized to the brane (the Casimir effect), which vanishes in the limit of 
large proper radius for the compact dimension. As will be discussed in a subsequent section the Casimir effect is finite 
and regulator independent in the limit that the cutoff is taken to infinity. Then by dimensional considerations if the 
fields on the brane are massless the right hand side must be of order [2-7rs(6)]~ ^ since the proper size of the compact 
dimension is the only scale in the problem. Then the equation above becomes 

|F(C3, b, a 2 ) - ^F(-c 2 , b, a 2 ) = const^ ^ l \\ - c 2 e-^ b ) 3 (l + c 2 e~i a - h r 5 ] , (34) 

where the constant above is of order one in units of the fundamental scale. While we see that both sides of this 
equation are the same order of magnitude even for large b we also see that b only appears in the combination e~z a2 . 
It is this combination which is determined in terms of the various tensions. Hence although the brane spacing is fixed, 
a large value for a 2 b is only possible through a fine tuning of parameters which is over and above the fine tuning 
necessary to set the four dimensional cosmological constant to zero. In the next section we shall show that by adding 
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a scalar field in the bulk a la Goldberger and Wise tt the matching conditions at the branes can be made more 
sensitive to b and the hierarchy can be made natural^]. 

Equating the various components of the metric across the boundary fixes the values of the coefficients /02, S02, 
/3, and S3. Once again there are four equations but there are now four unknowns so there is no further fine tuning 
required. 




Fig 1. Plot of warp factor of compact dimension (circumference of circles) versus r (vertical axis), illustrating why we 
call this the 'space needle' metric. The apparent singularity at the top, where r = 0, is a coordinate singularity only. 

This completes the determination of the metric in the absence of any bulk matter. The warp factor for the fifth 
dimension is plotted in figure 1 for a choice of parameters for which there is a large hierarchy. The warp factor for the 
usual 3+1 dimensions is qualitatively similar, except near r = where it goes to a constant. We call this the "space 
needle metric" for reasons which are obvious from the picture. 



B. Model with Bulk Scalar Fields 



In this section we show that when the model of the previous section is modified by the inclusion of bulk scalar 
fields, the hierarchy can be made natural. 
The action for the scalar field 

S M = J d 6 xV^G^(-d M i>d M i) - mV) + \f^GF(i))5{r - a) + y/^dH(ip)5(r - b) . (35) 

Here m is the mass of the scalar field in the bulk. We will be interested in the limit m 2 <§; a 2 , since this is where 
we naturally obtain a large hierarchy. The scalar field sources F(ip) and H(ip) are in general arbitrary functions. For 
simplicity we will take 



3 Note that the necessary fine tuning of the Planck 4-brane tension might also be natural if the 4-brane is an approximate 
BPS state of a nearly supersymmetric theory. Thus supersymmetry of the bulk action might provide an alternative solution to 
the hierarchy problem. 



G 



Fty) = Xxip (36) 
H(i>) = X 2 ^ • (37) 

The coupled equations for the gravity-matter system are difficult to solve. We shall address the problem by a 
successive approximation method. We will first solve for the gravitational field in the absence of the scalar field, 
assuming that the contribution to the stress tensor of this field is small. We will then solve for the scalar field in this 
geometry, and compute the stress tensor of the scalar field. We then can substitute this solution back into the 
Einstein equations to determine the correction to the geometry induced by the stress tensor of the scalar fielcj^]. If we 
wish, this procedure, which is essentially an expansion in T^/Af* 6 , can be carried out to higher orders but the lowest 
order will be sufficient to to fix the brane spacing, which we expect will receive only small corrections at higher orders. 

For simplicity we will take as a starting metric 

fo = 1 r < a (38) 
fo = e~ aa e ar a<r<b (39) 
f = e- aa e 2ab e- ar r>b , (40) 

s = r 2 r < a (41) 
s = a 2 e- aa e ar a< r < b (42) 
s = a 2 e- aa e 2ab e~ ar r>b . (43) 

This corresponds to a solution for the special case with no bulk cosmological constant for r < a but the same 
cosmological constant everywhere outside, and c 2 = 0. This metric is flat in the neighborhood of the origin and 
pure AdS outside a. This is a special case, with a\ and c 2 set to zero, of the more general class of metrics we have 
considered in section IIIA. Non-zero c 2 will be treated perturbatively in what follows and we will also indicate how 
to include ot\ perturbatively. Hence the only loss of generality arising from this starting metric is that we can only 
extend our conclusions to the more general class of metrics of the previous section when c 2 and a.\ are sufficiently 
small that a perturbative approach is valid. 

We first solve for the scalar field in this background metric. The equation of motion is 

- V" - Uj + ~ jj 4>' + = XiS(r - a) + X 2 6{r - b) . (44) 
The solutions in the three regions consistent with smoothness of tp at the origin and vanishing of if> at infinity are 

%j} = Ax (l + irnV + . . \ r < a (45) 

ip = A 2 e air + B 2 e a2r a<r <b (46) 
i/j = A 3 e CT3r r > b , (47) 

where we have neglected higher order terms in m 2 r 2 in equation (45), and 

(48) 
(49) 

V \ ^ / 

For positive m 2 a 2 is positive and a\ negative. Then 

03 = -er 2 ■ (50) 




4 We neglect quantum corrections to the stress tensor from the scalar field, as we do not expect these to qualitatively alter our 
conclusions. 
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The coefficients Ai and Bi are to be determined by matching the solutions for %p across the boundaries. We require 
continuity of ip across the boundaries and the following jump conditions 

A^'(o) = -F'{4>) (51) 

at the first boundary and 

AV/(6) = -H'ty) (52) 

at the second boundary. 

Since the expressions for the A^s and Sj's are complicated we will neglect effects of order e~ ab which are very small 
and further assume to 2 a 2 <C 1 so that such effects can also be neglected. These approximations will not affect our 
conclusions. Then in this limit the expressions for the Aj's and Bi's are 

A 2 = t — ~~~~r — e (^i)a e -a 2 b _ hl e -*i* (53) 

(C3 — (72)01 (71 

B 2 = -—**—e- nb (54) 

(0-3 - cr 2 ) 

A 3 = -—^—e-^ . (55) 

\(73 — (72 ) 

We now find the contribution to the stress tensor from the field tp, T$. 

16ttT/ = < -\{^' 2 + m 2 ) + F{iP)5(r - a) + H(^)5(r - b) > (56) 
167rV r = <i(^' 2 -m 2 )> . (57) 

Next, we substitute the stress tensor back into the Einstein equations to determine the corrections to the geometry. 
For convenience we define T = 87rT^,/2M» 4 , A = A/2M* 4 . The Einstein equations in the region a < r < b take the 
form 

The equations inside r < a can be obtained by setting a = in the equations above. The equations for r > b are 
identical to those above. We are interested in the correction to the geometry to linear order in TjM 2 . To obtain 
this we expand 

/ = /o(l + e) (61) 
s = «o(1 + k), (62) 

and linearize in e and k. Then for r < a the equations (58) and (60) above become 



2 

h" ~ 1*" - - 77- £ ' = ~T% ■ (64) 
2 2 r 2r 



e' = JrT6 6 (63) 



These equations yield 



e' = \rT\ (65) 

K '= l rT\-~ fdrr 2 T\, (66) 
* * r Jo 
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which can be integrated to get e and k. Here we have used the constraints of smoothness at the origin and the 
linearized form of the Bianchi identity. The latter is shown below. 

(rf\)' = Th . (67) 

For a < r < b the equations (59) and (60) when linearized become 

2e" + 5ae' = T\ (68) 
4ae' + cxk' = T% . (69) 



These yield 



e' = De-i ar + — - (70) 
5a 

k' = -4De~i ar + ^ , (71) 
5a 

which can be integrated to yield e and k. Here D is a constant of integration which must be determined from matching 
and once again the linearized form of the Bianchi identity for this region, shown below, has been used. 



In identical fashion we can get for the region r > b 



e' = *' = -^- (73) 

Here the requirement that the metric die away at infinity has been imposed. This completes the determination of 
bulk corrections to the metric. The final step is to match across the branes thereby determining the brane positions 
and the undetermined constant D. 

The conditions on the continuity of the metric across the brane are straightforward to satisfy because of the 
additional constants of integration that will be obtained on integrating the expressions for e' and k' . Once / has been 
normalized to one on the visible brane and the requirement of smoothness has been met at r = all these additional 
constants will have been fixed. 

We now focus our attention on the jump conditions on the derivatives. At the inner brane we have 

- ^Ae'(a) - lA«'(a) - \ (a - - \cx = - l -T^{a) (74) 

2Ae'(a) + 2a = - 7l 2 + ^A^(a) • (75) 

At the outer brane 

- |Ae'(6) - 1a«'(6) + \{2a) + \{2a) = /3 2 2 - ^(b) (76) 

2Ae'(6) - 4a = - 72 2 + ^j(b) . (77) 

These are four independent equations for the three unknowns a, b and D. Just as in the previous section, the 
metric and brane locations are completely determined with one overall fine tuning needed to find a Poincare invariant 
solution. 

First consider the situation without the scalar field, i.e. T = 0. This is the same case which was considered in the 
previous section. In general, as we found, the metric beween the branes need not be pure AdS even without a scalar 
field. The tensions of the branes must be finetuned to set D to zero. We are allowing these brane tensions (which are 
related to (3 2 and 7 2 in the above equations) to deviate from those fine tuned values by small amounts. Here small 
merely means that perturbation theory is valid. Then it is straightforward to verify that the equations (74) to (77) 
are simply linearized versions of equations (29) to (32) with a,\ — and the parameter D is proportional to c 2 . Hence 
even in the absence of the scalar field the metric between the branes is not of the AdS form, and our conclusions once 
the scalar field is introduced will be valid for this more general class of metrics. Setting T in the above formulas to a 
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constant non zero value for r < a corresponds to allowing non zero a\. Although it is straightforward to accomodate 
non zero ol\ perturbatively in this framework, for simplicity we shall keep a\ = in what follows below, and T will 
be related only to and will not include any piece from a cosmological constant. 

We now return to the general case with a scalar field. We are interested in a solution with all parameters of order 
one in units of the fundamental scale but we allow m 2 /a 2 to be somewhat less than one in order to obtain a hierarchy 

Rather than give the exact solution we will give the more informative order of magnitude results. For dimensional 
considerations any parameter of order one in units of the fundamental scale will be denoted by the appropriate power 
of a. 

The matching conditions at the inner brane, eqs. (74) and (75), determine a and D as functions of b. Provided 
m 2 b/a is less than or of order one, a and D are of order one in units of the fundamental scale. The dependence of D 
on b is 

D = 0(a)+0(m 2 b) . (78) 

This result will turn out to be crucial for the solution of the hierarchy problem. 

Both the equations at the outer brane are sensitive to b and a only through exponentially small terms. This is 
because the value of ip and its derivatives is exponentially insensitive to b and a in this region, which manifests itself 
in the forms of e' and n' close to the outer brane. This is similar to the insensitivity to b of equations (31) and (32). 
One might therefore worry that just as in the previous case an exponentially precise fine tuning will be necessary 
to get a hierarchy. However the finetuning in the previous case was related to the fact that near the Planck-brane 
the metric was very nearly pure AdS, which is a homogenous space, and the compact dimension was exponentially 
large, so that the location of the Planck brane had very little effect on the matching conditions at either brane. In 
the present case with a light bulk scalar the value of the scalar at the TeV brane depends more sensitively on b. To 
determine 6, add equations (76) and (77) to obtain 

b -De-i ab = (f3 2 f - ( 72 ) 2 • (79) 

This is the analogue of equation (33) in the previous section. Recall that (P2) 2 — arises from the vacuum 

energy of quantum fields localised on the brane and is of order ae~^ ah . Then this equation, together with equation 
((78|) which relates D to 6, determines b as 

6 = (4) . (80) 

Clearly m need not be much smaller than a to get a sizable hierarchy. The difference of equations (76) and (77), 
which we have not yet used, becomes the condition that the effective four dimensional cosmological constant is zero, 
which is the usual finetuning^]. 

One limitation of the above approach is that the inner brane necessarily has negative tension along the non compact 
directions. This is a consequence of the choice of a starting metric with zero a\ and C2 and the fact that deviations 
from this metric can only be perturbative. The nonzero ol\ case is not simple but there is a limit in which it is 
tractable — that in which (a\a) 2 is perturbatively less than one, even though a.\ is not small. In this limit, neglecting 
higher order terms in (air) we find 

f j = \a^r+... (81) 

S - = 2 --\ a 2 r + ... (82) 
s r 6 

We can substitute for this in equation (44) and find that to the order shown the solution (45) is unchanged. Now 
if a 2 a > a.2 the inner brane can have positive tension and it is straightforward to verify that all the other conclusions 
above go through as before. 



5 Work is in progress to see whether this finetuning can be made more natural though supersymmetry of the bulk action when 
the Planck brane is approximately BPS. 
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C. Physical Implications 



To determine the effective four dimensional Planck scale we concentrate on the higher dimensional Einstein action. 
When the two extra dimensions are integrated out this will contain the usual four dimensional Einstein action. 



S G = 2Mt J d 6 x M V^GR 6 ■ (83) 
Expanding g^ v = r]^ v + and integrating out the two extra dimensions we see that 

S G = 2M$ J d e x M y/^R 4 f^ + ... = Ml j d A x»^R A + ... , (84) 



and consequently 



2ttM a J drj^fs = M 2 . (85) 

Rather than do this integral exactly we will be satisfied with an estimate. The dominant contribution to the integral 
comes from the region close to the outer brane where the integrand is very large and the forms of / and s are very 
nearly simple exponentials. For the purposes of the estimate we set 013 = a 2 = a -1 = M* = a. 

Ml = 0(a 2 e? ab ) . (86) 

This large exponent is responsible for the hierarchy between the Planck scale and the weak scale. 

Next we analyze the spectrum of linearized tensor fluctuations to see if it is consistent with the results of gravitational 
experiments. We neglect the effect of the scalar field in what follows since we do not expect it to qualitatively affect 
our results because its contribution to the energy density is small. 

Expanding 

G>„ = fv^ + V ( 87 ) 
and substituting this into the Lagrangian we get the following equation for the fluctuation h in the bulk. 

-h"-~h'- (J^ h+^a 2 h=jm 2 h, (88) 

where m 2 = —r\ ilv V ll V v an d we are restricting our attention to modes that have no momentum in the compact 
direction. These are expected to be separated by a mass gap from the heavier modes with nonzero momentum in the 
fifth direction. 

The boundary conditions that h has to satisfy are h'(0) = at the origin and 

2A^ = - 7 2 . (89) 

Clearly h = f is a solution of this equation with m = by comparison with eqns. (5) and (6). This is the massless 
graviton. There is a continuum of other solutions for all positive to 2 , as can be seen from the asymptotic behavior of 
the equation. However to extract these solutions is not easy, because of the complicated forms of / and s. However 
since we are only interested in order of magnitude estimates we can simplify the problem by considering instead a 
simpler problem that retains the physics we are interested in. Consider the metric 

/ = 1 r <a (90) 
/ = e- aa e ar a<r<b (91) 
/ = e- aa e 2ab e- ar r > b (92) 

s = r 2 r <a (93) 
s = a 2 e- aa e ar a<r<b (94) 
s = a 2 e- aa e 2ab e- ar r>b . (95) 
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This corresponds to a solution for the case with no cosmological constant for r < a but the same cosmological 
constant everywhere outside. This metric is flat in the neighborhood of the origin, AdS outside a. It clearly has a 
form very similar to the metric we are interested in for r 3> a, although unlike our case, the positions of the branes 
are not fixed and one has negative tension. In what follows we will assume the main features of the excitations we 
are interested in are common to both metrics and proceed. 

With this approximation the equation in the region r < a where a is zero becomes 



_ h" - -ti = m 2 h . 



The solution of this equation to leading order in m 2 a is 



r 

2„2 



h = N 



-I 1 2 2 

1 m r 



(96) 



(97) 



where N is a normalization constant. Since we arc primarily interested in the light modes this will suffice. In the 
region between the branes the equation for the modes has the form 



2 2 / 

This equation admits a solution in terms of Bessel functions. The solution is 

h = Ne-^ ar {A 2 J i [mq 2 ] + B 2 J_s[mq 2 ]) 
where A 2 and B 2 are constants and q 2 is defined by 

2 e~i ar 

Q2 = 1 — • 

a e -2 aa 

The closed form expressions for the relevant Bessel functions are 



Ja(x) = \ — [ sin x 



J_s.(x) = \ — cos X 
" irx 1 



4-i 

X 1 

4-1 

X Z 



+ 3- 



^cosrr 

X 

sinx 



x 



(98) 



(99) 



(100) 



(101) 

(102) 
(103) 



Since we are interested in values of m such that mqi < 1 wc can conveniently approximate the solution between 
the branes by 



h = N 



A 2 — 



, 2m Y 
2 1 ~) 



1 



+ ^B 2 \ — 



2m 



Coming to the region r > b we can also obtain a solution in terms of Bessel functions. 

h = Ne+i ar (A 3 Js [mq 3 ] + B 3 J_,[mq 3 }^ 

where q 3 is defined by 



(104) 



(105) 



93 



a e -2 



aa+ab 



(106) 



The values of the A's and i3's are to be determined by matching. The boundary conditions to be satisfied are 
continuity of h across the various boundaries and the following jump conditions on the derivatives at r — a and r = b 
respectively 



A 



~h 

4 



-2a . 



(107) 
(108) 
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The constant N is to be determined by normalization. Since this is not quite in the form of a eigenvalue problem 
we make some simple transformations which render it so. 
Defining 



g = const h 



dr 



1 

77 



(109) 
(110) 



we get an equation for g as a function of q which has the form of an eigenvalue equation for m 2 with unit density 
function. For the continuum modes it is the outer region which is relevant for normalization. But here s and / are 
proportional so we can conveniently choose the constant so that g = h in this region. Also q and q 3 as defined differ 
at most by an additive constant. Hence for the continuum modes we merely have to normalize the solution for h for 

r > b with respect to q 3 . 

Matching at the inner brane we find that A 2 is of order yjm/a, and B 2 is of order (m/a)?. Then matching at the 
outer brane we find that 



A 3 = 
B 3 = 



say 
(")' 



In the far asymptotic region the A 3 mode is dominant. Normalizing to a box of size L we find 

1 m 



N = 



—= — e 

VL a 



The situation is slightly different for the zero mode. We write the normalizable wave function as 

Now the integral relevant for normalization 

J dqgl = N% J dr-Lfhi = N 2 J drfsi . 



(Ill) 
(112) 

(113) 
(114) 
(115) 



This is the same integral that appears in determination of the four dimensional Planck scale. By exactly the same 
methods we obtain, on normalizing to unity 



N Q = 0(ae~* ab ) 



(116) 



We are now in a position to determine the corrections to gravity from the Kaluza Klein excitations of the graviton. 
The change in the potential energy between two masses mi and mi on our brane is given by 



AV = 



I 



dm 



N 2 
Nfa 



Le 



O 



Gm\mi- 



e -|a6 



O 



10 



-32 



r 3 (TeV) 3 



(117) 



From this it is clear that deviations from Newtonian gravity are highly suppressed at long distances. 

We now explain why we expect this model to have the same physical implications as the model we started out with. 
Essentially for r 3> a the general solution of both models will have similar form. The only difference will be in the 
magnitudes of the coefficients Ai and B<i. Although these coefficients are determined by matching in the interior their 
order of magnitude follows from simple dimensional considerations. This then implies that A 3 and B 3 and hence the 
normalization of the modes can be fixed by dimensional considerations. Hence the two theories will give the same 
order of magnitude estimate for the the corrections to Newtonian gravity. 

Because of the isometry of the compact dimension this theory contains a massless "gravi-photon" — a Kaluza-Klcin 
U(l) gauge boson. However no light or massless fields will carry non-trivial U(l) charge since they have no momentum 
in the compact dimension. We expect that other than the modes we have already discussed, the remaining spectrum 
of gravitational excitations will be massive. 

A more comprehensive study of the phenomenological implications of this model is left for future work. 
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D. Stress Energy Tensor for a Field Localized to a Brane 



In this section we consider the form of the stress energy tensor for a field localized to a brane having the metric 
Gab = diag(—l, 1, 1, 1, 1) but in which the fifth dimension is compact and has proper size a. We show that in the 
ground state the stress tensor does indeed have the form given in eqn. (23). For simplicity we limit ourselves to the 
case of a free scalar field. 

The Lagrangian has the form 

L = -^d A 4>d A 4>-^m 2 ^ 2 . (118) 

We are interested in the expectation value of the stress tensor in the ground state. 

< Tf >=< - G B A L > . (119) 

Because the ground state possesses translational symmetry 

< Tf >= I J d'x < T? > (120) 

where V is the volume of the four space dimensions. Performing a Fourier expansion for the field and making use 
of the canonical commutation relations this reduces to 

< T o°> =EEiVk 2 + fc5^ + m 2 (121) 

k k 5 

<k> = -EE ^7 TPrfpT^ ■ (122) 

Since the usual three space dimensions are infinite the sums over momenta in these three directions can be replaced 
by integrals. However since the fifth dimension cp is compact the momenta in this direction remain discrete, k§ = n/a 
where n is an integer. 



< T o° > = E / d 3 Wk 2 + hk 5 + m 2 



(123) 



<g>-g^/ rt - vTOTg - (124) 

These integrals are infinite and must be regulated in order to yield sensible physical results. We will use a Pauli- 
Villars regulator, adding massive fields with appropriate statistics until all the divergences have been removed. (We 
could get similar results in a theory with spontaneously broken supersymmetry.) We simplify to the special case where 
the boson field is massless. Then all the divergences can be removed by adding three fields with opposite statistics 
having masses M,M and 2M and two fields with the same statistics which both have mass \/3M. Here M is assumed 
to be some kind of cut off for the theory. 

Performing the integrals and adding the contributions from the various fields we get the finite but regulator depen- 
dent results 

T m = 5 ;E i^)^^^^) + 2(kl + 3M 2 ) 2 ln(kl + 3M 2 ) - 2(fc 5 2 + M 2 ) 2 ln(k 2 5 + M 2 ) - (k 2 5 + 4M 2 ) 2 ln(k 2 5 + AM 2 )] (125) 

T " = - E (^^MPs) + 2(A: 5 2 + ZM 2 )ln{kl + 3M 2 ) - 2(k 2 5 + M 2 )ln(k 2 5 + M 2 ) - (fc 5 2 + AM 2 )ln{k 2 b + AM 2 )] (126) 

Tl = . (127) 

All other components of the stress energy tensor vanish, and it clearly has the form of eqn (23). In a supersymmetric 
theory the scale M will be related to the scale of supersymmetry breaking. We now estimate and Xf in various 
limits. 

When M <C (1/a) we can approximate and T| to leading order in (1/a) 2 as 0(M /a) and 0(M 2 /a 3 ) respectively. 
Clearly in this limit 1 and 7 2 are not equal. 
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The relative difference between Tf and T" must be finite when the cutoff M is taken to infinity and must vanish 
as the size of the fifth dimension becomes infinite. We now estimate this difference as a function of a and M when 
the dimensionless quantity aM is large. 

To do this we attempt to replace the sum we are interested in by a sum of integrals. We begin by observing that 
the integral below can be broken up into a sum of integrals over equal subdomains. 

dk 5 T(k 5 ) =J2 dp 5 T(p 5 ) . (128) 

fc 5 Jk s 

Here T represents an arbitrary function of k$. If the function T is smooth it can be Taylor expanded 

T(p 5 ) = T(k 5 ) + T'(k 5 )(p 5 - k 5 ) + T"(fc 5 ) (p5 ~ fc5)2 + . . . (129) 



Then performing the integrals over the subdomains we get 



f dfc 5 T(fc 5 )=]T 

J fe5 



lnh) + ^nk 5) + ±_T»(k 5) + 



(130) 



Now if T is T,™ or Tf the first term on the right has the form we are interested in. Also notice that the other 
terms on the right then involve fewer powers of k§ in the numerator and hence for the seventh term and beyond the 
sums for individual fields are finite and straightforward to estimate. 

T can also be chosen to be derivatives to arbitrary order of TJ£ or Tf . Since these quantities ocur in the expansions 
for T£ or Tf they can then be substituted back to obtain systematic expansions for and Tf in terms of integrals. 

A complication that arises for the case of a massless field is that the fourth derivative and higher of T£ and Tf are 
not well defined at k§ — 0. We account for this by separating this point from the sum and approximating the rest of 
the sum by integrals from (1/a) to infinity and — (1/a) to n egat ive infinity. 



To obtain a reasonable estimate we must expand in eq. (129) to at least seventh order in order to account for all 
possible divergences as powers or logarithms of M. The calculation is straightforward but lengthy and the details will 
not be presented here. The result is that the difference between T^ and Tf is finite and of order (1/a) 5 in the limit 
that the cutoff M is large. 

This (1/a) 5 result for the difference could have been anticipated. Since the only counterterm allowed by general 
covariance is the cosmological constant which contributes equally to both and Tf , the difference between these 
two must be finite and regulator independent in the limit that the cutoff is taken to infinity. For a massless field, a 
is the only available dimensionful parameter. This result is just a higher dimensional form of the Casimir force. 



IV. CONCLUSIONS 



We have constructed a set of solutions to Einstein's equations in six or more dimensions, and exhibited a six 
dimensional set up, "the space needle" , with two concentric positive tension 4-branes, which each have one compact 
dimension. Gravity is mostly localized to the outer brane while we assume the standard model lives on the inner 
brane, explaining the apparent weakness of gravity in our world. There are no massless moduli associated with either 
the size of the compact dimension or the brane locations. This provides an explicit demonstration that the gauge 
hierarchy problem can be solved in six dimensions, without supersymmetry, and with negligible corrections to gravity 
at distances longer than an inverse TeV. Gravitational effects do become strong at energies of order a TeV. We leave 
discussion of gravitational collider phenomenology of new noncompact dimensions for future work. 

We do not address the important issue of how to obtain chiral fermions on the standard model brane. Ordinary 
dimensional reduction by compactifying the fifth dimension on a circle always results in a non-chiral theory. One 
simple alternate possibility is to have the standard model live on a 3-brane at the center of space. Then it is only 
necessary to have one 4-brane — the Planck brane. The metric is simply the a — > limit of the space needle metric. 
Alternatively, it may be possible to generalize our mechanism to additional dimensions with some compactification 
which does allow a chiral effective theory on the TeV brane. 

We also do not address the cosmological constant problem. The effective four dimensional cosmological constant 
depends on a complicated function of the bulk and brane parameters, and may be finetuned to zero or to a small 
acceptable value. 
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